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ABSTRACT 

One of the challenges in hyperspectral data analysis is the 
presence of mixed pixels. Mixed pixels are the result of low 
spatial resolution of hyperspectral sensors. Spectral unmix¬ 
ing methods decompose a mixed pixel into a set of endmem- 
bers and abundance fractions. Due to nonnegativity constraint 
on abundance fraction values, NMF based methods are well 
suited to this problem. In this paper multilayer NMF has been 
used to improve the results of NMF methods for spectral un¬ 
mixing of hyperspectral data under the linear mixing frame¬ 
work. Sparseness constraint on both spectral signatures and 
abundance fractions matrices are used in this paper. Evalu¬ 
ation of the proposed algorithm is done using synthetic and 
real datasets in terms of spectral angle and abundance angle 
distances. Results show that the proposed algorithm outper¬ 
forms other previously proposed methods. 

Index Terms — Hyperspectral data, spectral unmixing, 
multilayer NMF, sparseness constraint. 

1. INTRODUCTION 

Hyperspectral sensors are an effective tool for remote sensing 
applications. These sensors offer high spectral resolution, so 
they can gather a lot of information about the spectral content 
of materials that are present in the scene. But the main draw¬ 
back of these sensors is the low spatial resolution that cause 
mixed pixels to appear in hyperspectral images |[T]. The con¬ 
cept of mixed pixels in hyperspectral images is illustrated in 
Fig. [TJusing a toy example. 

There are three different categories of spectral unmixing 
methods: geometrical, statistical and sparse regression meth¬ 
ods Q. Vertex component analysis (VCA) Q is one of the 
classic methods in the geometrical category that is used in this 
paper as a baseline of the proposed method. 

Another class of algorithms that are used for spectral un¬ 
mixing purposes are methods based on nonnegative matrix 
factorization (NMF). NMF method gets a lot of attention in 
this area of research because of nonnegativy constraint on 
abundance fraction values. Different constraints can be used 
in NMF methods, for example in Q minimum volume con¬ 
straint is used or in Q graph regularized constraint has been 
used. In this paper multilayer NMF Q is used to unmix 



Fig. I: Toy example of mixed pixels concept in hyperspectral 
images Q. 

hyperspectral data. In this approach the observation matrix 
will be decomposed in different layers. So the final estimated 
spectral signatures can be modeled as a product of sparse ma¬ 
trices obtained in the layers. The reason why decomposing 
a general matrix into sparse matrices will improve the results 
of NMF method is still an open problem ||^ . Experiments on 
synthetic and real datasets show that the proposed method can 
achieve better results in comparison with other methods. 

The rest of the paper are as follows. Section [^introduces 
mathematical description of linear mixing model for solving 
spectral unmixing problem. Section [^briefly reviews NMF 
method and commonly used sparseness constraint on abun¬ 
dance fractions matrix. In section [^ multilayer structured 
NMF has been described to use in hyperspectral unmixing. 
Experiments on synthetic and real datasets has been done to 
evaluate the proposed method in section]^ Einally section]^ 
concludes the paper. 

2. LINEAR MIXING MODEL (LMM) 

Spectral unmixing problem can be solved under linear or non¬ 
linear frameworks In this paper linear mixing model 
(LMM) has been used. The mathematical representation of 
this model is expressed in Q. 

X = AS + N, (1) 

where X G refers to observation matrix, L denotes 

number of spectral bands and N denotes total number of pix¬ 
els. A S and S G refer to signature and 


















abundance fractions matrices respectively. P denotes num¬ 
ber of pure signatures (called endmembers) in the scene and 
N S refers to measurement noise. Two constraints on 

abundance fraction values should be considered: Abundance 
Nonnegativity Constraint (ANC) and Abundance Sum to one 
Constraint (ASC). These constraints are expressed in © and 

0 . 

> 0 , ( 2 ) 

P 

^S,, = l. (3) 

i=l 

3. NONNEGATIVE MATRIX EACTORIZATION 
(NMF) 


Considering the sparseness constraint on spectral sig¬ 
natures and abundance fraction matrices, the proposed cost 
function for each layer is expressed in 

Omlnmf = ||X - ASIII + a||A||i/2 + A||S||i/2, (8) 

where a and A are regularization parameters and control the 
impact of sparseness constraints on spectral signatures and 
abundance fractions respectively. 

Optimizing the cost function in using multiplicative 
rules results in estimated A and S. Note that the cost func¬ 
tion in (j^ should be minimized in each level of the algorithm. 
In this paper to consider the ASC constraint on abundance 
fraction values, FCLS method GD has been used. 


Nonnegative matrix factorization can decompose a matrix 
into a product of two matrices efficiently. NMF can be ex¬ 
pressed mathematicaly as: 

X « AS. (4) 


5. EVALUATION RESULTS 

To evaluate the performance of the proposed algorithm, spec¬ 
tral angle distance (SAD) and abundance angle distance 
(AAD) are used. These metrics are defined in (|^ and ([TOll. 


To calculate this approximation, the following cost func¬ 
tion can be dehned using Euclidean distance fTO) : 


SADrui = cos 




(9) 


= ||X-AS||| (5) 

As explained in section[T] due to ANC on abundance frac¬ 
tions matrix, NMF methods has been widely used for spectral 
unmixing of hyperspectral data. Without employing other 
constraints to NMF method, the results are not satisfactory. 
One of the most common constraints is sparseness constraint 
on abundance fraction matrix GD- Considering this con¬ 
straint the cost function of NMF method changes to the fol¬ 
lowing equation: 

OsparseNMF = |1X — AS|||. -f A||S|ji/ 2 , (6) 

4. MULTILAYER NMF 

Cichocki and Zdunek |j7) proposed multilayer structure for 
NMF to reduce the risk of getting stuck in local minima and 
improve the performance of NMF methods. In this paper we 
have used this idea to solve spectral unmixing of hyperspec¬ 
tral data. 

Using multilayer structure, the first layer of algorithm is 
a basic decomposition of observation matrix into Ai and Si 
matrices. In the second layer Si will be decomposed into A 2 
and S 2 and this process will be repeated to reach the maxi¬ 
mum number of layers. Then the estimated spectral signature 
matrix can be calculated using Q. 

A = A 1 A 2 ... (7) 

where L^ax is the maximum number of layers. 


AAD a, =cos II ) (10) 

SAD and AAD are dehned for one endmember and one 
pixel respectively. To have an overall measure for all end- 
members and pixels, root mean square values of these metrics 
have been used in this paper and dehned in GU and G3- 

1 P 

rmsSAD = (11) 

i=l 

r ms AAD = ( — ^(AADaJ^)'/' (12) 

i=l 

5.1. Synthetic Dataset 

In this experiment, synthetic data has been generated using 
spectral signatures of USGS spectral library JT^. Selected 
materials from this library are shown in Fig. HThe process 
described in Q has been used in this paper to generate simu¬ 
lated data. 

Fig. shows the results of applying proposed method on 
synthetic dataset. This hgure illustrates the results in terms 
of rmsSAD and rmsAAD in comparison with VGA Q and 
L 1 / 2 -NMF pT) for different values of SNR ranging from 
IbdB to 45(167 As it can be seen from results the proposed 
method can unmix data more effectively compared to other 
methods. 





Fig. 2: Selected Materials from USGS Library. 



Fig. 3: Comparison of Results in terms of rmsSAD (Top) 
and rmsAAD (Bottom) for VCA, L 1 / 2 -NMF and MLNMF 
against SNR. 

5.2. Real Dataset 

In the second experiment Cuprite Nevada dataset collected by 
AVIRIS sensor 03 ) is used as real dataset for evaluating the 
proposed method. This dataset contains 250 x 191 pixels and 
188 spectral bands. This area consists of different minerals. 


Table 1: Comparison between methods in terms of rmsSAD 
for Cuprite Nevada dataset 


Method 

VCA SISAL L 1 / 2 -NMF MLNMF 

rmsSAD 

0.1047 0.1237 0.1138 0.0981 



Fig. 4; Estimated signatures (dotted in red) in comparison 
with USGS library signatures (continous in blue) for Cuprite 
dataset 


Reference spectra to compare the results are selected from 
USGS library (jT^. 

The proposed method has been applied to this dataset and 
results are compared with VCA Q , SISAL p3] and L 1 / 2 - 
NMF flT) in terms of rmsSAD andpresented iiiTabIe[T] Also 
in Fig.^0 and spectral signatures and abundance fraction 
maps that are extracted using the proposed method are illus¬ 
trated. Results conhrms that this method can effectively un¬ 
mix hyperspectral data. 


6. CONCLUSION 

Hyperspectral data contains mixed pixel due to low spatial 
resolution of the sensors. Spectral unmixing tries to decom¬ 
pose a mixed pixel into the spectral signature and abundance 
fractions. NMF methods has been widely used to solve spec¬ 
tral unmixing of hyperspectral data. In this paper to improve 
the performance of NMF method multilayer structure is used. 
Sparsity constraints on spectral signatures and abundance 
fractions have been added to the cost function. Synthetic 
and real datasets are used to evaluate the performance of the 
proposed mathod in comparison with other methods. Results 
in terms of SAD and AAD shows that the proposed method 
can effectively unmix the hyperspectral data. 
































































































Buddingtonite GDS85 D-206 


Alunite GDS82 Na82 


Andradite WS487 Garnet 


Kaolincmect H89-FR-5 


Kaolincmect KLF508 


Chalcedony CU91-6A 


Dumortierite HS190.3B 


Montmorillonite+Illite CM37 


Muscovite IL107 


Sphene HS189.3B 


Nontronite NG-l.a 


Andradite WS474 


Fig. 5 : Estimated abundance fraction maps for Cuprite dataset 
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